Anomalous c-axis charge dynamics in copper oxide materials 
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Within the t-J modeh the c-axis charge dynamics of the copper oxide materials in the underdoped 
and optimally doped regimes is studied by considering the incoherent interlayer hopping. It is shown 
that the c-axis charge dynamics is mainly governed by the scattering from the in-plane fluctuations. 
In the optimally doped regime, the c-axis resistivity is a linear in temperatures, and shows the 
metallic-like behavior for all temperatures, while the c-axis resistivity in the underdoped regime 
is characterized by a crossover from the high temperature metallic-like to the low temperature 
semiconducting-like behavior, which are consistent with experiments and numerical simulations. 
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The copper oxide materials are among the most com- 
plex systems studied in condensed matter physics. The 
complications arise mainly from (1) strong anisotropy in 
the properties parallel and perpendicular to the Cu02 
planes which arc the key structural element in the whole 
copper oxide superconducting materials, and (2) extreme 
sensitivity of the properties to the compositions (stoi- 
chiometry) which control the carrier density in the Cu02 
plane^. After over ten years of intense experimental 
study of the copper oxide materials, a significant body 
of reliable and reproducible data has been accumulated 
by using many probes, which indicates that the normal- 
state properties in the underdoped and optimally doped 
regimes are quite unusual in many aspects suggesting 
the unconventional metallic state realized^. Among the 
striking features of the normal-state, the quantity which 
most evidently displays the anisotropic property in the 
copper oxide materials is the charge dynamics^, which is 
manifested by the optical conductivity and resistivity. It 
has been shown from the experiments that the in-plane 
charge dynamics is rather universal within the whole cop- 
per oxide materials^'"^. The in-plane optical conductiv- 
ity for the same doping is nearly materials independent 
both in the magnitude and energy dependence, and shows 
the non-Drude behavior at low energies and anomalous 
midinfrared band in the charge-transfer gap, while the 
in-plane resistivity pab(T) exhibits a linear behavior in 
the temperature in the optimally doped regime and a 
nearly temperature linear dependence with deviations at 
low temperatures in the underdoped regime^"''. By con- 
trast, the magnitude of the c-axis charge dynamics in the 
underdoped and optimally doped regimes is strongly ma- 
terials dependent, i.e., it is dependent on the species of 
the building blocks in between the Cu02 planes^. Al- 
though the c-axis charge dynamics is very complicated, 
some qualitative features seem to be common, such as 
(1) for the optimally doped YBa2Cu3 06-i-a: and over- 
doped La2-2:Sr^Cu04 systems^'^, the transferred weight 
in the c-axis conductivity decaies as —>■ 1/lu at low en- 



ergies, which is in accordance with the metallic-like c- 
axis resistivity Pc{T) for all temperatures, and (2) for 
the lower doping the temperature dependent c-axis re- 
sistivity Pc{T) is characterized by a crossover from the 
high temperature metallic-like behavior to the low tem- 
perature semiconducting-like behavior^"^. The nature of 
the c-axis charge dynamics in the copper oxide materials 
is of great importance, as the superconducting mecha- 
nism is closely associated with the anisotropic normal- 
state properties^. 

Since the undoped copper oxide materials are antifer- 
romagnetic Mott insulators, upon hole doping, the anti- 
ferromagnetic long-range-order is rapidly destroyed and 
the unusual metallic state emerges^*^. In this case, many 
researchers believe that the essential physics is contained 
in the doped antiferromagnets"'^^'^^, which may be effec- 
tively described by the t-J model acting on the space 
with no doubly occupied sites^"'^^^^. On the other hand, 
there is a lot of evidence from the experiments and nu- 
merical simulations in favour of the t-J model as the ba- 
sic underlying microscopic model^*^'^^. Within the two- 
dimensional (2D) t-J model, the in-plane charge dynamic 
of the copper oxide materials has been extensively stud- 
ied theoretically as well as numerically-'^''' ^^. Since the 
understanding of the charge dynamics in the copper ox- 
ide materials is not complete without an understanding of 
the c-axis charge dynamics, therefore a number of alter- 
native mechanisms for the c-axis charge dynamics have 
been proposed^^'^^, and the most reliable results for the 
c-axis charge dynamics have been obtained by the numer- 
ical simulation'^^. To shed light on this issue, we, in this 
paper, apply the fermion-spin approach'^®'^^ to study the 
c-axis charge dynamics based on the t-J model by consid- 
ering the interlayer coupling. Within each Cu02 plane, 
the essential physics properties are described by the 2D 
t-J model as, 
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where f] — ±aox, iaoj/, ^o is the lattice constant of 
the square planar lattice, which is set as the unit here- 
after, i refers to planar sites within the 1-th Cu02 plane, 
C/jg. (C/icr) are the electron creation (annihilation) op- 
erators, S/i — Cj^aCii/2 are the spin operators with 
a ~ [gx, cfy,az) as the Pauli matrices, and ^ is the chem- 
ical potential. The Hamiltonian (1) is supplemented by 
the on-site local constraint, ^^Cl^^Cua < 1, i.e., there 
be no doubly occupied sites. For discussing the c-axis 
charge dynamics, the hopping between Cu02 planes is 
considered as^* 



H^-tcJ2 CLC'+'7c*^ + h.c. + J2 Hh 
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where fjc — ±coz, cq is the interlayer distance, and has 
been determined from the experiments^'^ as cq > 2ao. In 
the undcrdopcd and optimally doped regimes, the exper- 
imental results^^^ show that the ratio R = pc{T) / pab{T) 
ranges from R ^ 100 to i? > 10^, this large magnitude 
of the resistivity anisotropy reflect that the c-axis mean 
free path is shorter than the interlayer distance, and the 
carriers are tightly confined to the Cu02 planes, and also 
is the evidence of the incoherent charge dynamics in the 
c-axis direction. Therefore the c-axis momentum can not 
be defined^^. Moreover, the absence of the coherent c- 
axis charge dynamics is a consequence of the weak in- 
terlayer hopping matrix element tc, but also of a strong 
intralayer scattering, i.e., tc <C i, and therefore the com- 
mon Cu02 planes in the copper oxide materials clearly 
dominate the most normal-state properties. In this case, 
the most relevant for the study of the c-axis charge dy- 
namics is the results on the in-plane conductivity (Tab (^) 
and related single-particle spectral function A(k,uj). 

Based on the 2D t-J model, the self-consistent mean- 
field theory in the undcrdopcd and optimally doped 
regimes has been developed^" within the fermion-spin 
approach^^, which has been applied to study the pho- 
toemission, electron dispersion and electron density of 
states in the copper oxide materials, and the results are 
qualitative consistent with the experiments and numeri- 
cal simulations. Moreover, the in-plane charge dynamics 
in the copper oxide materials has been discussed^^ by 
considering the fluctuations around this mean-fleld solu- 
tion, and the results exhibits a behavior similar to that 
seen in the experiments* and numerical simulations^^. 
In the fermion-spin theory^^'^*^, the constrained electron 
operators in the t-J model is decomposed as. 
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with the spinless fermion operator hi keeps track of the 
charge (holon), while the pseudospin operator Si keeps 
track of the spin (spinon) , then it naturally incorporates 
the physics of the charge-spin separation. The main ad- 
vantage of this approach is that the electron on-site local 



constraint can be treated exactly in analytical calcula- 
tions. In this case, the low-energy behavior of the t- 
J model (2) in the fermion-spin representation can be 
rewritten as""^**, 

H = tcY. hU^MSt^S^+^^, + S^S+^^J + Y,Hi, (4a) 
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where Jeff ~ J[(l — (5)^ — 
parameter (j) = {huhu+f)), 



6^], the holon particle-hole 
and 5;^ and S^^ are the 



pseudospin raising and lowering operators, respectively. 
These pseudospin operators obey the Pauli algebra, i.e., 
they behave as fermions on the same site, and as bosons 
on different sites. It is shown^'' that the constrained elec- 
tron operator in the t-J model can be mapped exactly 
onto the fermion-spin transformation defined with an ad- 
ditional projection operator. However, this projection 
operator is cumbersome to handle for the actual calcu- 
lation possible in 2D, we have dropped it in Eq. (4). It 
has been shown in Ref.^^ that such treatment leads to 
the errors of the order 5 in counting the number of spin 
states, which is negligible for small doping 5. 

In the framework of the charge-spin separation, an 
electron is represented by the product of a holon and a 
spinon, then the external field can only be coupled to one 
of them. loffe and Larkin^'^ and Li et al.'^'^ have shown 
that the physical conductivity (j{ijj) is given by. 
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where a^^^{Lu) and (T*^*'(cj) arc the contributions to 
the conductivity from holons and spinous, respectively. 
Within the Hamiltonian (4), the c-axis current densi- 
ties of holons and spinous are given by the time deriva- 
tive of the polarization operator using Heisenberg's equa- 
tion of motion as, j^ == '^tcexJ^ifj^iVchj^fi^^hu and 

jc'^ = tce(t)c T,lf|,^ f]c{S^,S-_^f^^^ + S-^S+_f^J, respectively, 
where tc — tcXc/x is the effective interlayer holon hop- 
ping matrix element, and the mean-field spinon and 
holon order parameters are defined^" as Xc = i^H^i+fi i)' 
X = {SuS^t+f;)^ and (/)c = (huhi+fiat)- As in the pre- 
vious discussions^*, a formal calculation for the spinon 
part shows that there is no the direct contribution to 
the current-current correlation from spinous, but the 
strongly correlation between holons and spinous is con- 
sidered through the spinon's order parameters entering 
in the holon part of the contribution to the current- 
current correlation, therefore the charge dynamics in the 
copper oxide materials is mainly caused by the charged 
holons within the Cu02 planes, which are strongly renor- 
malized because of the strong interactions with fluctua- 
tions of the surrounding spinon excitations. In this case, 



the c-axis optical conductivity is expressed^^ as (Tc{iLj) — 
— Imllc {llj)/u! with the c-axis holon current-current cor- 
relation function n[^\t - t') = {{jc'\t)ji''\t'))). In the 
case of the incoherent charge dynamics in the c-axis di- 
rection, i.e., the weak interlayer hopping tc ^ t, this 
c-axis holon current-current correlation function He (cu) 
can be evaluated in terms of the in-plane holon Green's 
function g{k,uj)^'^'^^, then we obtain the c-axis optical 
conductivity as^'*'^^^^^, 



o'c(w) = ■^i^icexco)'^ 



1 ^^ f°° dcu' .^ , 
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vanishes at higher temperature (T > 0.4 J). These re- 
sults are also qualitative consistent with the experimental 
results^'^ and numerical simulations^*. 

The quantity which is closely related with the c-axis 
conductivity is the c-axis resistivity pdT), and can be 
expressed as, 
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This c-axis resistivity has been evaluated numerically and 
the results at the doping S = 0.12 and 5 = 0.06 for the 
parameters t/J = 2.5, tc/t = 0.04, and cq/qq = 2.5 are 
shown in Fig. 3(a) and Fig. 3(b), respectively. In the 
underdoped regime, the behavior of the temperature de- 



where npii-o) is the Fermi distribution functions, the in- 
plane holon spectral function Ah{k,u!) = — 2Im(7(fc,a;), 
while the in-plane holon Green's function g(k,uj) has 
been obtained by considering the second-order correction 
for holons due to the antiferromagnetic fluctuations, and 
given in Ref.^^. As pointed in Ref.^*, the approximation 
assumption of the independent electron propagation in 
each layer has been used in the above discussions, and is 
justified for tc <C t, therefore the c-axis conductivity is 
essentially determined by the properties of the in-plane 
spectral function. We have performed a numerical cal- 
culation for the c-axis optical conductivity crc(ct;), and 
the results at the doping 5 = 0.12 (solid hne), 5 — 0.09 
(dashed line), and 5 = 0.06 (dot-dashed hne) for the pa- 
rameters t/J = 2.5, ic/t = 0.04, and co/ao = 2.5 at 
the temperature T=Q are shown in Fig. 1, where the 
charge e has been set as the unit. From Fig. 1, it is 
found that crc(a;) is composed of two bands separated at 
Lo ~ 0.4t, the higher-energy band, corresponding to the 
"midinfrared band" in the in-plane optical conductivity 
fah(w)^'^*, shows a broad peak at a; ~ 0.7t, moreover, 
the weight of this band is strongly doping dependent, 
and decreasing rapidly with dopings, but the peak posi- 
tion does not appreciably shift to higher energies, which 
is consistent with the experimental results^^^. On the 
other hand, the transferred weight of the lower-energy 
band forms a sharp peak at w < 0.4t, which can be 
described formally by the non-Drude formula, and our 
analysis indicates that this peak decay is — > 1/lo at low 
energies as in the case of crQf,(w)'*'^^. In comparison with 
o'ah('^)^'*, the present results also show that the values 
of (Tc(w) are by 2 ~ 3 orders of magnitude smaller than 
those of Uabi'^) in the corresponding energy range. For 
further understanding the property of (Tc(a;), we have also 
discussed the finite temperature behavior of i7c(w), and 
the numerical results at the doping 5 — 0.12 for the pa- 
rameters t/J = 2.5, tc/t = 0.04, and co/ao — 2.5 with 
T = 0.2J (solid hne) and T = 0.5J (dashed line) are 
plotted in Fig. 2, which show that adio) is temperature 
dependent for w < 1.2i and almost temperature inde- 
pendent for LO > 1.2t, while the higher-energy band is 
severely suppressed with increasing temperatures, and 
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FIG. 1. The c-axis optical conductivity at the doping 

5 = 0.12 (solid line), 5 = 0.09 (dashed line), and 5 = 0.06 
(dot-dashed line) for t/J = 2.5, ic/t = 0.04, and co/ao = 2.5 
with the temperature T — 0. 
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FIG. 2. The c-axis optical conductivity at the doping 
S = 0.12 for t/J = 2.5, ic/t = 0.04, and Co/ao = 2.5 with 
the temperature T — 0.2J (solid line) and T = 0.5J (dashed 
line). 



pendence of Pc{T) shows a crossover from the high tem- 
perature metaUic-like {dpc{T)/dT > 0) to the low tem- 
perature semiconducting-Uke {dpc{T)/dT < 0), but the 
metalhc-hke temperature dependence dominates over a 
wide temperature range. In comparison with the in- 
plane resistivity pab{TY'^, it is shown that the crossover 
to the semiconducting- like range in Pc{T) is obviously 
linked with the crossover from the temperature linear to 
the nonlinear range in pab{T), and are caused by the 
pseudogap observed in the normal-state, but Pab{T) is 
only shghtly affected by the pseudogap^^, while Pc{T) 
is more sensitive to the underlying mechanism. Our re- 
sults also show that there is the common origin for these 
crossovers. Therefore in this case, there is a general trend 
that the copper oxide materials show nonmetallic Pc{T) 
in the underdoped regime at low temperatures. While in 
the optimally doped regime, Pc{T) is a linear in tempera- 
tures, and shows the metallic-like behavior for all temper- 
atures. These results are qualitative consistent with the 
experimental results^"^ and numerical simulation^^. It 
has been shown from the experiments^^ that the charge 
dynamics in some strongly correlated ladder materials 
shows the similar behaviors. 

In the above discussions, the central concerns of the 
c-axis charge dynamics in the copper oxide materials are 
the two dimensionality of the electron state and incoher- 
ent hopping between the Cu02 planes, and therefore the 
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FIG. 3. The c-axis electron resistivity at the parameter 

t/J = 2.5, tc/t — 0.04, and co/ao ~ 2.5 for (a) the doping 
6 = 0.12 and (b) <5 = 0.06. 



c-axis charge dynamics in the present fermion-spin pic- 
ture is determined by the in-plane charged holon fluctua- 
tions. In this case, the c-axis scattering rate is associated 
with the in-plane scattering rate, and can be roughly de- 
scribed by the imaginary part of the self-energy of the 
charged holons within the Cu02 planes, which is consis- 
tent with the "dynamical dephasing" theory proposed by 
Leggett^^. 

In the fermion-spin theory ^^, the charge and spin de- 
grees of freedom of the physical electron are separated as 
the holon and spinon, respectively. Although both holons 
and spinous contributed to the charge and spin dynam- 
ics, but it has been shown that the scattering of spinous 
dominates the spin dynamics^^, while the results of the 
in-plane charge dynamics^^ and present c-axis charge dy- 
namics shows that scattering of holons dominates the 
charge dynamics, therefore the notion of the charge-spin 
separation naturally accounts for all the qualitative fea- 
tures of the normal-state properties of the copper oxide 
materials. To our present understanding, the main rea- 
sons why the fermion-spin theory based on the charge- 
spin separation is successful in studying the normal-state 
property of the strongly correlated copper oxide mate- 
rials are that (1) the electron single occupancy on-site 
local constraint is exactly satisfied in the analytic calcu- 
lation. Since the anomalous normal-state property of the 
copper oxide materials are caused by the strong electron 
correlation in these systems^°'^^, and can be effectively 
described by the t-J model^°~^^, but the strong electron 
correlation in the t-J model manifests itself by the elec- 
tron single occupancy on-site local constraint, then the 
satisfaction of this local constraint is equivalent to that 
the strong electron-electron interaction has been properly 
treated. This is why the crucial requirement is to treat 
this constraint exactly in the t-J model in the analytic 
discussions. (2) Since the local constraint is satisfied even 
in the mean-field approximation within the fermion-spin 
theory ^^, the extra gauge degree of freedom related to 
the common " flux" phase problem occurring in the slave- 
particle approach^^ does not appear here, which is con- 
firmed by our previous discussions within the mean-field 
theory^", where the photoemission, electron dispersion 
and electron density of states in the copper oxide materi- 
als have been studied, and the results are qualitative sim- 
ilar to that seen in the experiments and numerical simu- 
lations. (3) As mentioned above, the dropping the pro- 
jection operator in Eq. (4) will only lead to errors of the 
order S in counting the number of spin states within the 
common decoupling approximation^^. This because that 
the constrained electron operators Cia in the t-J model 
can be also mapped onto the slave-fermion formulism^^ as 

Cia — hjaia with the local constraint ^j^i + X^o- ^L^*'^ ~ 
1. We can solve this constraint by rewritting the bo- 
son operators aia in terms of the CP^ boson operators 
b^a as Qicr = (1 - hlh^Y / "^ bi„ « (1 - h\h^/2)bi^ supple- 
mented by the local constraint X^ct ^L^io- — 1- Since 
the CP^ boson operators 6i| and hii with the local con- 



straint can be identified with the pseudospin lowering 
and raising operators in the fermion-spin approach^®, re- 
spectively, then the spinon propagator in the restricted 
Hilbert space can be written as Dii{i — j,t — t') = 
{{[1 ~ hl{t)h,{t)/2]; [I - h\(t')K{t')/2]))D{^ '-j,t- t') « 
[1 - (5 - 0{6^)]D{i - j,t- i'), where D{i - j,t - t') is 
the spinon propagator within the fermion-spin approach. 
In this case, the extra spin degrees of freedom in the 
fermion-spin theory only lead to the errors of the order 
5 in calculating the spinon propagator within the com- 
mon decoupling approximation^^, which is negligible for 
small doping 5. This is why the theoretical results of the 
spin dynamics within the fermion-spin approach^* are 
qualitative consistent with the experiments and numeri- 
cal simulations. 

In summary, we have studied the c-axis charge dynam- 
ics of the copper oxide materials within the t-J model by 
considering the incoherent interlayer hopping. Our re- 
sults show that the c-axis charge dynamics is mainly gov- 
erned by the scattering from the in-plane charged holon 
fluctuations. The c-axis optical conductivity and resistiv- 
ity have been discussed, and the results are qualitative 
consistent with the experiments and numerical simula- 
tions. 

Finally we also note that since the structure of the 
building blocks in between the Cu02 planes for the 
chain copper oxide materials is different from these for 
the nochain copper oxide materials, some subtle differ- 
ences between the chain and nochain copper oxide ma- 
terials for the c-axis charge dynamics have been found 
from the experiments^^®. It has been suggested* that 
for the nochain copper oxide materials the doped holes 
may introduce the disorder in between the Cu02 planes, 
contrary to the case of the chain copper oxide materi- 
als, where the increasing doping reduces the disorder in 
between the Cu02 planes. It is possible that the dis- 
order introduced by the doped holes residing between 
layers in the nochain copper oxide materials in the un- 
derdoped regime may modify the interlayer hopping ele- 
ments, which leads to the subtle differences between the 
chain and nochain copper oxide materials for the c-axis 
charge dynamics. These and other related issues are un- 
der investigation now. 
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